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Abstract 

A model is developed for a detailed investigation of the current flowing 
through a cylindrical nanosize MOSFET with a close gate electrode. The 
quantum mechanical features of the lateral charge transport are described by 
Wigner distribution function which is explicitly dealing with electron scatter- 
ing due to acoustic phonons and acceptor impurities. A numerical simulation 
is carried out to obtain a set of I-V characteristics for various channel lengths. 
It is demonstrated that inclusion of the collision term in the numerical simula- 
tion is important for low values of the source-drain voltage. The calculations 
have further shown that the scattering leads to an increase of the electron 
density in the channel thereby smoothing out the threshold kink in I-V char- 
acteristics. An analysis of the electron phase-space distribution shows that 
scattering does not prevent electrons from flowing through the channel as a 
narrow stream, and that features of both ballistic and diffusive transport may 
be observed simultaneously. 

I. INTRODUCTION 

During the last decade significant progress has been achieved in the scaling of the metal- 
oxide-semiconductor field-effect transistor (MOSFET) down to semiconductor devices with 
nanometer sizes. In Ref . |^ , the fabrication of silicon quantum wires with lengths and widths 
of about 60 nm and ~ 20 nm respectively is reported. The conductance of those quantum 
wires was measured for a wide range of temperatures, from 25 to 160 K. The fabrication and 
the investigation of a 40 nm gate length n-MOSFET are reported in Ref. |^. The resulting 
nanosize n-MOSFET operates rather normally at room temperature. Using nano imprint 
lithography, a field effect transistor (FET) with a 100 nm wire channel was fabricated 
and the characteristics of this FET at room temperature were investigated. However, the 
small size of nanoscale MOSFET with a wide Si substrate negatively influences the device 
characteristics due to the floating body effect. Short channel effects together with random 
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effects in the silicon substrate are very well known to cause a degradation of the threshold 
voltage and the appearance of uncontrollable charge and current in regions far from the gate 
electrode. Therefore in a nanoscale conventional MOSFET, the controlling ability of the 
gate electrode is substantially weakened. Recently, considerable attention has been paid to 
SOI (Si-on-insulator) MOSFETs, which are prospective for creating new nanosize devices. 
In Refs. @-|^, a MOSFET with very thin SOI was theoretically investigated on the basis 
of a 2D analytical model, while in Ref. |^, a ID model was used. Drain- induced barrier 
lowering was considered and the physical mechanisms which determine the subthreshold 
slope (S-factor) were analyzed [^. As a result, a substantial reduction of the short-channel 
effect in the SOI MOSFET as compared to that in the bulk devices was established. As 
was shown theoretically in Ref. [§, the use of a lightly doped source and drain leads to 
an increase of the effective channel length what allows one to weaken the drain-induced 
barrier lowering. In an SOI MOSFET with a Si-Ge source 0, an improved drain-to-source 
breakdown voltage is achieved due to the absorption of excess holes in the channel region. In 
a transistor device with a channel sandwiched between oxide layers (dual-gate MOSFET), 
the floating body effects are significantly suppressed [jTO[ . A theoretical model of a dual-gate 
device is described in Ref. [0. In Ref. [0] we have investigated the thermal equilibrium 
state of a nanoscale cylindrical silicon-based MOSFET device with a close gate electrode 
(MOSFETCGE). An advantage of the latter is the complete suppression of the floating body 
effect caused by external influences. Moreover, the short-channel effect in these devices can 
be even weaker than that in a dual-gate structure. 

The main goal of the present work is the investigation of quantum transport in a nanosize 
MOSFETCGE device. We have developed a flexible 2D model which optimally combines 
analytical and numerical methods and describes the main features of the MOSFETCGE 
device. The theoretical modeling of the quantum transport features involves the use of the 
Wigner distribution function formalism |ll3|-|l5|. The paper is organized as follows. In Sec- 
tion II, a description of the system is presented in terms of a one-electron Hamiltonian. In 
Section III the quantum Liouville equation satisfied by the electron density matrix is trans- 
formed into a set of one dimensional equations for partial Wigner distribution functions. A 
one-dimensional collision term is derived in Section IV. In Section V we describe a numerical 
model to solve the equations which have been derived for the partial Wigner distribution 
function. In Section VI the results of the numerical calculations are discussed. Finally, in 
Section VII we give a summary of our results and conclusions about the influence of the 
scattering processes in nanosize MOSFETs. 



II. THE HAMILTONIAN OF THE SYSTEM 

We consider a cylindrical nanosize MOSFETCGE structure (Fig. 1) described by cylin- 
drical coordinates {r,(j),z), where the z-axis is chosen to be the symmetry axis. In the 
semiconductor pillar the electron motion is determined by the following Hamiltonian 

where ^(r) = ^(1") + ^(1") is the potential energy associated with the energy barrier and 
the electrostatic field, respectively; mj- and m^! are the effective masses of the transverse (in 
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(x, ?/)-plane) and longitudinal (along z-axis) motion of an electron of the j'-tli valley. 
The electrostatic potential energy satisfies Poisson's equation 

AV;(r) = — (-n(r) + Noir) - iV^(r)) , ^ = 1, 2, (2) 

where Si and 62 are the dielectric constants of the semiconductor and oxide layers, re- 
spectively; n(r), A''£)(r)and A^a(i') are the concentrations of electrons, donors and acceptors 
respectively. In our calculations we assume that the source electrode is grounded whereas 
the potentials at the drain and gate electrodes are equal Vds and Vq, respectively. 

The study of the charge distribution in the cylindrical nanosize MOSFETCGE structure 
in the state of the thermodynamical equilibrium (see Ref. [0) has shown that the concen- 
tration of holes is much lower than that of electrons so that electron transport is found to 
provide the main contribution to the current flowing through the MOSFET. For that reason, 
holes are neglected in the present transport calculations. 



III. THE LIOUVILLE EQUATION 

In this section we consider ballistic transport of electrons. Neglecting scattering processes 
and inter-valley transitions in the conduction band, the one-electron density matrix can be 
written as 

P(r,r')=EP.(r,r'), (3) 

j 

where pj (r, r') is the density matrix of electrons residing in the j-th valley satisfying Liou- 
ville's equation 

th^ = [H„p,]. (4) 

In order to impose reasonable boundary conditions for the density matrix in the electrodes, 
it is convenient to describe the quantum transport along the z-axis in a phase-space rep- 
resentation. In particular, we rewrite Eq. (^) in terms of ( = {z + z')/2 and rj = z — z' 
coordinates and express the density matrix pj as 

r+oo 

p,(r,r)= J2 TT dke''''^f,msm's'iC,k)^jmsir±,z)^*^,^,{r'^,z'), (5) 

ms,m's' 



with a complete set of orthonormal functions \Efjms(r±, -z)- According to the cylindrical 
symmetry of the system, these functions take the following form: 

^,m.(r±, z) = -^^^^,{r, z)e'"'t (6) 



The functions ipjms{r, z) are chosen to satisfy the equation 
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2mj 



I d f d \ 
r dr \ dr I 



ipjms{r, z) + F(r, z)ipj^s{r, z) = £jms{z)iJjms{r, z), (7) 
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which describes the radial motion of an electron. Here Sjms{z) are the eigenvalues of Eq. 
for a given value of the 2;-coordinate which appears as a parameter. It will be shown, that 
£jms{z) plays the role of an effective potential in the channel, and that \l/jms(r±, -z) is the 
corresponding wavefunction of the transverse motion at fixed z. Substituting the expansion 
(I) into Eq. (|^), and using Eq. (|^), we arrive at an equation for fjmsm's'iC: k) : 

'KT fjmsm's'iCj k) ~^ ^ [ W^jmsm's' (C) ^ k ) fjmsm' s' {C y ^ )dk 

"J 

+00 



dt ji^C'^'"""'*'^' ' h 

3 —00 
+00 

^ j M;;1,,,(C, K k')f,ms,m's[{C, k')dk', (8) 



where 



+00 

,AC, k-k') = — I (£:,„,(C + v/2) - 8,^,s'{C, - v/2)) e'^^'-'^^drj, (9) 
zm J 



+00 



M;::L'AC k, k') = ^ J [^.^^f „,,,(C + v/^, y) + ^..^f^,,, (C - r//2, k')] e^^^'-'^^^dv, (10) 



^mssi^^jk^ 11 ^jmssi(^) ~l~ ||Cjmssi(^) I ^ q/- ~I~ 2/u j , (H) 



II 

2m ji 




and 



brmss^iz) = I ^IJ*^,{r, z)—'^jmsi{r, z)rdr, (12) 



92 
dz 

/d 

Note that Eq. (|^) is similar to the Liouville equation for the Wigner distribution function, 
which is derived to model quantum transport in tunneling diodes (see Ref. 0]). The first 



drift term in the right-hand side of Eq. (P) is derived from the kinetic-energy operator of 
the longitudinal motion. It is exactly the same as the corresponding term of the Boltz- 
mann equation. The second component plays the same role as the force term does in the 
Boltzmann equation. The last term in the right-hand side of Eq. (|^) contains the operator 

^jmsm's'iCy k, k'), which mixes the functions fjmsm's' with different indexes s,s'. It appears 
because ipjmsir, z) are not eigenfunctions of the Hamiltonian ([l|). The physical meaning of 
the operator Mj^^^,g,{(, fc, k') will be discussed below. 

In order to solve Eq. (^, we need to specify boundary conditions for the functions 
fjmsm's'iC^k). For a weak current, electrons incoming from both the source and the drain 
electrodes, are assumed to be maintained in thermal equilibrium. Comparing Eq. (^ with 
the corresponding expansion of the density matrix in the equilibrium state, one obtains the 
following boundary conditions 
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fjmsm's'iO, k) = S 

ss'^mm' 

2 [exp {EjsmkP - EfsP) + 1]" , k>0, 

fjmsm's'{L,k) = 6ss'Smm'2[exp{EjsmkP - EfdP) + 1]'^ , k < 0, (14) 



where the total energy is Ejsmk = fi^k'^/2mj + Sjsm{0) for an electron entering from the 
source electrode {k > 0) and Ejsmk = h'^k'^/2m}j + Sjsm{L) for an electron entering from 
the drain electrode {k < 0). Eps and Epp, are the Fermi energy levels in the source and 
in the drain, respectively. Note, that Eq. (|T^ meets the requirement of imposing only one 
boundary condition on the function f jmsm's' {(, k) at a fixed value of k as Eq. (|D is a first 
order differential equation with respect to C- Generally speaking, the solution of Eq. (H) 



with the conditions (|T^) depends on the distance between the boundary position and the 
active device region. Let us estimate how far the boundary must be from the active device 
region in order to avoid this dependence. It is easy to show that the density matrix of 
the equilibrium state is a decaying function oi t] = z — z' . The decay length is of the 

order of the coherence length At = at high temperature and of the inverse Fermi 



wavenumber k^ = /2m}jEp at low temperature. So, it is obvious, that the distance 
between the boundary and the channel must exceed the coherence length or the inverse 
Fermi wavenumber, i. e. L ^ \p or L ^ kp^. For example, at T = 300 K the coherence 
length At ~ 3 nm is much less than the source or drain lengths. 

The functions f jmsm's' {C^^)^ which are introduced in Eq. (|^), are used in calculations 
of the current and the electron density. The expression for the electron density follows 
directly from the density matrix as n(r) = p(r,r). In terms of the functions fjmsm's'iCy k), 
the electron density can be written as follows: 

^ +00 

^('^) = ^ J2 J f jmsm's' iz,k)dk<iIjmsir±,z)<iJ*j^,^,{r±,z). (15) 

jmsm's'_^ 



It is well-known [|l^], that the current density can be expressed in terms of the density matrix 

(16) 



yr-^ eh f d d \ , 

j r,t = ^ : — - — r,r',t 

^ 2m j% \d\: or' j 



The total current, which flows through the cross-section of the structure at a point z, can 
be obtained by an integration over the transverse coordinates. Substituting the expansion 
(^) into Eq. ([1^) and integrating over r and <y9, we find 



J — e - J dk ^fjmsmsi^^k) IP Cjmss'i^) J dk liafj^gj-ns' {z, k) , (17) 



hk 2eTi 

' jmsms\Z ) ) IP 

j,m,s ^" rilj nij j.m 

s' 

where Im/ is the imaginary part of /. The first term in the right-hand side of Eq. (|l] 
is similar to the expression for a current of the classical theory [jT6[. The second term, 
which depends on the non-diagonal functions f jmsms' only, takes into account the effects of 
intermixing between different states of the transverse motion. 

The last term in the right-hand side of Eq. (^ takes into consideration the variation of the 
wavefunctions ipjms{f, z) along the z-axis. In the source and drain regions, the electrostatic 
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potential is essentially constant due to the high density of electrons. In these parts of the 
structure, the wavefunctions of the transverse motion are very weakly dependent on z, and 
consequently, the operator Mj^^j^,^, has negligible effect. Inside the channel, electrons are 
strongly localized at the Si/Si02 interface as the positive gate voltage is applied. Earher 
calculations, which we made for the case of equilibrium , have shown that in the channel 



the dependence of 'ipjms{f,z) on z is weak, too. Therefore, in the channel the effect of the 

operator M^mj^'s' negligible. In the intermediate regions (the source-channel and the 
drain-channel), an increase of the contribution of the third term in the right-hand side of 
Eq. (H) is expected due to a sharp variation of ipjms{f,z). Since Mj^sm's' couples functions 
fjmsm's'iCyk) with different quantum numbers (jms), it can be interpreted as a collision 
operator, which describes transitions of electrons between different quantum states of the 
transverse motion. Thus, the third term in the right-hand side of Eq. (H) is significant only 
in the close vicinity of the p-n junctions. Therefore, this term is assumed to give a small 
contribution to the charge and current densities. Under the above assumption, we have 
treated the last term in the right-hand side of Eq. (steady state of the system in a zeroth 
order approximation with respect to the operator M^l^sm's'- Neglecting the latter, one finds 
that, due to the boundary conditions (|14|), all non-diagonal functions fjmsm's'iC: ^) ("^ 



or s 7^ s') need to be zero. 

In the channel, the energy of the transverse motion can be approximately written in the 
form [r^l 



where Sjsiz) is the energy associated with the radial size quantization and h'^m'^ /2mj-gR*l is 
the energy of the angular motion with averaged radius R*g. Hence, in Eq. ^ for the diagonal 
functions fjmsmsiC^ the difference SjmsiC + ^/2) — £jms{.C ^ ''7/2) can be substituted by 
£js{C + ^7/2) — £js{C, — ^7/2). Furthermore, summation over m in Eq. (^) gives 

/,.(C, ^) - ^ y WjsiC, k - fc')/,s(C, k')dk' = (19) 



rrij d( 



with 



fjsiCik) — fjmsmsjCi k)- (20) 



In Eq. ([T9|) the following notation is used 



+00 



WjsiC, k) = -^ I {£,s{C + - ^i.(C - W2)) sm{kr^)d7^. (21) 

— CO 

The effective potential Sjs{z) can be interpreted as the bottom of the subband {j,s) in 
the channel. The function fjs{C-,k) is referred to as a partial Wigner distribution function 
describing electrons which are travelling through the channel in the inversion layer subband 
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IV. ELECTRON SCATTERING 



In this section we consider the electron scattering from phonons and impurities. For this 
purpose we introduce a Boltzmann-hke single collision term ||T^, which in the present case 
has the following form 



fjsmk ^ ^ i^PjsmkJ's'm'k' fj's'm'k' Pjsmk,j' s'm'k' f jsmk) 

j's'm'k' 



(22) 



As was noted above, we have neglected all transitions between quantum states with 
different sets of quantum numbers j and s. In the source and drain contacts the distribution 
of electrons over the quantum states of the angular motion corresponds to equilibrium. 
Consequently, due to the cylindrical symmetry of the system, we may fairly assume that 
across the whole structure the electron distribution is given by 



fjsmk{z) = fjs{z, k)w 



jsmi 



where 



U!jsm 



2™2 



Ph m 



(23) 



(24) 



J "IS, 



is the normalized Maxwellian distribution function with respect to the angular momentum 
m. The integration of the both sides of Eq. (^) over the angular momentum gives the 
one-dimensional collision term 



St fjs{z, k) — ^ [Pjs[k, k )fjg[z, k ) Pjsik , k)fjg[z, k)) , 
k' 



(25) 



where 



Pjsik, ) — ^ ^ Pjsmk,jsm'k''^jsm' ■ (2^) 
mm' 

This collision term is directly incorporated into the one-dimensional Liouville equation (|TU|) 



as 



W,s{z, k, k') = W,s{z, k-k') + Pjsiz, k, k') - 6k,k' E Pjs{z^ k', k), (27) 

k' 



where Wjs{z, k, k') is the modified force term in Eq. (plQl). 

In this work we consider scattering by acceptor impurities and acoustic phonons described 
by a deformation potential. The scattering rates are evaluated according to Fermi's golden 
rule 



jsmkjsm'k' 



27r 

T 



{jsm'k'l Hint \jsmk) S {Ejsm'k' - Ejsmk) , 



(28) 



where Hint is the Hamiltonian of the electron-phonon or the electron-impurity interaction. 
Hereafter, we model the potential of an ionized acceptor as U{r) = 4Tre'^Rl/ei6{r), where Rs 
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determines a cross-section for scattering by an impurity. Consequently, the absolute value 
of the matrix element is 



\{jsm'k'\U{r - r^) \jsmk)\ = Ane Rjei^J.^ri^Zi). 



(29) 



Averaging this over a uniform distribution of acceptors results in the following scattering 
rate 



R 



pi 

jsmk, jsm'k' 



Ci / ip%{r, z)6 {E jsm'k' - Ejsmkk) rdr, 



(30) 



where Ci = Na (Ane'^R'^/ei)'^ /Ti and Na is the acceptor concentration. 

At room temperature the rate of the scattering by acoustic phonons has the same form. 
Indeed, for T = 300 K the thermal energy fc^T ^ huj^^, therefore the acoustic deformation 
potential scattering is approximately elastic, and the emission and absorption rates are equal 
to each other. For low energies we can approximate the phonon number as N^^ ~ ksT/huq S> 
1 and the phonon frequency = voq, where vq is the sound velocity. Assuming equipartition 
of energy in the acoustic modes, the scattering rate is 



:>ph 



27r 



jsmkjsm'k' = Y^Ph \ {jsm'k'\ e*^-"" \jsmk) 5 {E jsm'k' 



E 



jsmk) 1 



(31) 



where the parameter Cph = AT^ksT /'^^Tipvlfi. Integrating over q yields the scattering rate 
,y in the form (|30|) with Cph instead of Ci. The full scattering rate Pjsmkjsm'k' = 

pt _|_ pPh 

jsmk, jsm'k' ' ' 

scattering rate 



pph 

jsmk, jsm 

jsmk,jsm'k' "T J- jsmk,jsm'k' ^i^Gu inserted into Eq. in order to obtain the one-dimensional 



where 



Pjsiz,k,k') = iCi + Cph)aiz)F 



h'k 



21,2 



2m; 



K'k 
2m}] 



(32) 



aiz 



2mi 



R 



R,s{z) / ij%{r, z)rdr, F{x) = e-''/'Ko{\x\/2), 



where Kq(x) is a McDonald function ||T^. In calculations of the scattering by acoustic 
phonons the following values of parameters for Si are used: S =9.2 eV, p = 2.3283-10'^ kg/m^, 
vo = 8.43 ■ 10^ cm/s 



V. NUMERICAL MODEL 



The system under consideration consists of regions with high (the source and drain) 
and low (the channel) concentrations of electrons. The corresponding electron distribution 
difference would produce a considerable inaccuracy if we would have attempted to directly 
construct a finite-difference analog of Eq. (0). It is worth mentioning that, in the quasi- 
classical limit, i. e. SjsiC + ^/2) — SjsiC ~ vf^) ~ ^^^^) Eq. (|TP|) leads to the Boltzmann 
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equation with an effective potential which has the following exact solution in the equilibrium 
state: 



f!siCk) = -Y: 



27.2 



exp 



2m] 



./5 + ^,,,(C)/? + 



EfI3\+1 



(33) 



For numerical calculations it is useful to write down the partial Wigner distribution function 
as fjs{(,k) = fjsiCyk) + ffsiCyk). Inserting this into Eq. (0), one obtains the following 
equation for //^(C, k): 



^§c^'U(' j ^^s{c. k - k')ff^{c, k')dk' = i?,,(c, k), 



where 



m] dC 



BUC, ^) = ^ / 7 dv UsiC + |) - ^..(C - |) - ^^v) sin [{k - k')v] fJ^iC, k'). 
—00 —00 ^ ' 

(35) 

The unknown function ffXC-, k) takes values of the same order throughout the whole system, 
and therefore is suitable for numerical computations. In the present work, we have used the 
finite-difference model, which is described in Ref. ||T3|. The position variable takes the set 
of discrete values Ci = for = 0,...,Ai'^}. The values of k are also restricted to 
the discrete set kp = {2p — A^^ — l)Ak/2 for {p = 1, . . . ,Nk}. On a discrete mesh, the 
first derivative ^^{Q,kp) is approximated by the left-hand difference for kp > and the 



(34) 



right-hand difference for kp < 0. It was shown In Ref. [|T^, that such a choice of the finite- 
difference representation for the derivatives leads to a stable discrete model. Projecting the 
equation (^) onto the finite-difference basis gives a matrix equation L-f = b. In the matrix 
L, only the diagonal blocks and one upper and one lower co-diagonal blocks are nonzero: 



(36) 





( ^1 


-E 


. 


. 


\ 




-V 


A2 


-E . 


. 




L = 





-V 


^3 . 


. 






[ 





. 




-J 



Here, the Nk x matrices Ai, E, and V are 

2m]AC 



[Ailppf — 5pp, 



h\2p - Nk - l)Ak 



^jsiCiy kp k 



P )i 



[E]pp, = 5ppS[^^-p], [V]pp, = 5pp,e[p-^^], (37) 



and the vectors are 
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[fi\p = fjsiCi, kp), and [bi]p = Bj^iQ, K)^ ^ = 1' ^C-i? ^ = 1' ^k- 



(38) 



A recursive algorithm is used to solve the matrix equation L • f = b. Invok- 
ing downward elimination, we are dealing with Bi = {Ai — V Bi_i)^^ E and Ni = 
{Ai — VBi_i)^^ {bi + VNi_^i) [i = 1, . . . , A''^) as relevant matrices and vectors. Then, up- 
ward elimination eventually yields the solution /j = -Bi/i+i + {i = — 1, . . . , 1). If an 
index of a matrix or a vector is smaller than 1 or larger than A*"^ — 1, the corresponding term 
is supposed to vanish. 

In the channel, the difference between effective potentials Sjs{C) with different (j, s) is of 
the order of or larger than the thermal energy fc^T. Therefore, in the channel only a few low- 
est inversion subbands must be taken into account. In the source and drain, however, many 
quantum states (j, s) of the radial motion are strongly populated by electrons. Therefore, 
we should account for all of them in order to calculate the charge distribution. Here, we can 
use the fact that, according to our approximation, the current flows only through the lowest 
subbands in the channel. Hence, only for these subbands the partial Wigner distribution 
function of electrons is non-equilibrium. In other subbands electrons are maintained in the 
state of equilibrium, even when a bias is applied. So, in Eq. ( P^ for the electron density, we 
can substitute functions fjmsms{z, k) of higher subbands by corresponding equilibrium func- 
tions. Formally, adding and subtracting the equilibrium functions for the lowest subbands 
in Eq. (|T5|), we arrive at the following equation for the electron density 



hoo 



n{r) = neq{r) + ^ H J dk fj,{z, k) \il'js{r, z)\^ - f-^{z, k) ij]%r, z) , (39) 

where neg(r) and ip'jl{r, z) are the electron density and the wavefunction of the radial motion 
in the state of equilibrium, respectively. The summation on the right-hand side of Eq. (RW) is 
performed only over the lowest subbands. Since the electrostatic potential does not penetrate 
into the source and drain, we suppose that the equilibrium electron density in these regions 
is well described by the Thomas-Fermi approximation: 

ne,(r) = Nc^Fy^ {(3 (e\/(r, z) + Ep - Ec)) , (40) 

where the Fermi integral is 

Here Nc is the effective density of states in the conduction band and Ep is the Fermi level 
of the system in the state of equilibrium. 



VI. NUMERICAL RESULTS 



During the device simulation three equations are solved self-consistently: (i) the equation 
for the wavefunction of the radial motion (^, (ii) the equation for the partial Wigner distri- 
bution function (^) and (iii) the Poisson equation (Q). The methods of numerical solution 
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of Eqs. (1^ and are the same as for the equihbrium state [|T2| . The numerical model for 



Eq. (|T9|) was described in the previous section. In the present calculations the four lowest 
subbands (j = 1,2 s = 1,2) are taken into account. The electron density in the channel 
is obtained from Eq. ([15D , whereas in the source and drain regions it is determined from 
Eq. (0). The calculations are performed for structures with a channel of radius R = 50 nm 
and for various values of the length: Let =40, 60, 70 and 80 nm. The width of the oxide 
layer is taken to be 4 nm. All calculations are carried out with A*"^ = 100 and A^^ = 100. The 
partial Wigner distribution functions, which are obtained as a result of the self-consistent 
procedure, are then used to calculate the current according to Eq. (|TCp. 

We investigate two cases: ballistic transport and quantum transport. The scattering 
of the electrons is taken into account. The distribution of the electrostatic potential is 
represented in Fig. 2 for Vds = 0.3 V and Vg = 1 V. This picture is typical for the MOSFET 
structure, which is considered here. The cross-sections of the electrostatic potentials for r = 
0, 30, 40, 45, 48, 50 nm are shown in Fig. 3. The main part of the applied gate voltage falls in 
the insulator (50 nm < r < 54 nm). Along the cylinder axis in the channel, the electrostatic 
potential barrier for the electron increases up to about 0.4 eV. Since the potential along the 
cylinder axis is always high, the current mainly flows in a thin layer near the semiconductor- 
oxide interface. This feature provides a way of controlling I^g through the gate voltage. 
Varying V^s and Vq mainly changes the shape of this narrow path, and, as a consequence, 
influences the form of the effective potential Sjs{z). As follows from Figs. 2 and 3, the 
radius of the pillar can be taken shorter without causing barrier degradation. At the p-n- 
junctions (source-channel and drain-channel) the electrons meet barriers across the whole 
semiconductor. These barriers are found to persist even for high values of the applied 
source-drain voltage and prevent an electron flood from the side of the strongly doped 
source. The pattern of the electrostatic potentials differs mainly near the semiconductor- 
oxide interface, where the inversion layer is formed. In Fig. 4 the effective potential for the 
lowest inversion subband (j = 1, s = 1) is plotted as a function of z for different applied bias 
Vds = 0, . . . , 0.5 V, Vg = 1 V, Lch = 60 nm. It is seen that the effective potential reproduces 
the distribution of the electrostatic potential near the semiconductor-oxide interface. In the 
case of ballistic transport (dashed curve), the applied drain-source voltage sharply drops 
near the drain-channel junction (Figs. 3 and 4). The scattering of electrons (solid curve) 
smoothes out the applied voltage, which is now varying linearly along the whole channel. 
Note, that the potential obtained by taking into account scattering is always higher than that 
of the ballistic case. The explanation is clear from Fig. 5, where the linear electron density 
is plotted for Vg = 1 V and L^h = 60 nm. It is seen that, due to scattering, the electron 
density in the channel (solid curve) rises and smoothes out. Hence, the applied gate voltage 
is screened more effectively, and as a result, the potential exceeds that of the ballistic case. 
It should be noted that at equilibrium (Vds = 0) the linear density and the effective potential 
for both cases (with and without scattering) are equal to each other. This result follows 
from the principle of detailed balance. 

The current-voltage characteristics (the current density / = J/27tR vs. the source- 
drain voltage Vds) are shown in Figs. 6 and 7 for the structures with channel lengths Lch = 
40, . . . , 80 nm. At a threshold voltage Vds ~ 0.2 V a kink in the I-V characteristics of 
the device is seen. At subthreshold voltages Vds < 0.2 V the derivative dVds/dJ gives the 
resistance of the structure. It is natural, that scattering enhances the resistance of the 



11 



structure (solid curve) compared to the ballistic transport (dashed curve). Scattering is also 
found to smear the kink in the I-V characteristic. At a voltage Vds > 0.2 V a saturation 
regime is reached. In this part of I — V characteristics, the current through the structure 
increases more slowly than it does at a subthreshold voltage. The slope of the I — V curve in 
the saturation regime rises when the length of the channel decreases. This effect is explained 
by a reduction of the p-n junction barrier potential as the length of the channel becomes 
shorter than the p-n junction width. In Fig. 7, one can see that, when the transistor is 
switched off {Vq < 0.5 V), the influence of the scattering on the current is weak. This fact 
is due to a low concentration of electrons, resulting in a low amplitude of the scattering 
processes. 

In Figs. 8a and 8b the contour plots of the partial Wigner distribution function (j = 
1, s = 1) are given for both cases (a - without and b - with scattering). The lighter regions 
in these plots indicate the higher density of electrons. Far from the p-n-junction, where 
the effective potential varies almost linearly, the partial Wigner distribution function can be 
interpreted as a distribution of electrons in the phase space. When electrons travel in the 
inversion layer without scattering, their velocity increases monotonously along the whole 
channel. Therefore, in the phase-space representation the distribution of ballistic electrons 
looks as a narrow stream in the channel (Fig. 8a). As it is expected, scattering washes out 
the electron jet in the channel (see Fig. 8b). It is worth mentioning that the electron stream 
in the channel does not disappear. It means that in this case the electron transport through 
the channel combines the features of both diffusive and ballistic motion. 

VII. SUMMARY 

We have developed a model for the detailed investigation of quantum transport in MOS- 
FET devices. The model employs the Wigner distribution function formalism allowing us to 
account for electron scattering by impurities and phonons. Numerical simulation of a cylin- 
drical nanosize MOSFET structure was performed. I — V characteristics for different values 
of the channel length were obtained. It is shown that the slope of the I — V characteristic 
in the saturation regime rises as the channel length increases. This is due to the decrease of 
the p-n junction barrier potential. 

Finally, we have demonstrated that the inclusion of a collision term in numerical sim- 
ulation is important for low source-drain voltages. The calculations have shown that the 
scattering leads to an increase of the electron density in the channel and smoothes out the 
applied voltage along the entire channel. The analysis of the electron phase-space distribu- 
tion in the channel has shown that, in spite of scattering, electrons are able to flow through 
the channel as a narrow stream although, to a certain extent, the scattering is seen to wash 
out this jet. Accordingly, features of both ballistic and diffusive transport are simultaneously 
encountered. 
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FIGURE CAPTIONS 



Fig. 1. Scheme of the cyhndrical nanosize MOSFET. 

Fig. 2. Distribution of the electrostatic potential in the MOSFET with Lch — 60 nm at 

VG = iy and Vds = 0.3 V. 

Fig. 3. Cross-sections of the electrostatic potentials without scattering (dashed curves) and 
with scattering from acceptor impurities and from an acoustic deformation potential (solid 
curves) . 

Fig. 4. Effective potential as a function of z for various Vds, Lch ~ 60 nm. 

Fig. 5. Linear electron density in the channel as a function of z for various Vds, Lch = 60 nm. 

Fig. 6. Current-voltage characteristics at Vg = 1 V for different channel lengths. 

Fig. 7. Current-voltage characteristics for MOSFET with Lch = 40 nm. 

Fig. 8. Contour plots of the partial Wigner distribution function fjs{z, k) for the lowest 

subband (j = 1, s = 1) at Vg = 1 V, Vds — 0.3 V, Lch — 60 nm: a - without scattering, b - 

with scattering. 
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